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Abstract 



We study a phenomenological model for weak ferromagnetic Fermi liquids and investigate the properties of 
^ ■ the spin waves in the model. The Landau kinetic equation is used to derive, in addition to the known Goldstone 

O ' mode, a new spin-wave mode - the first Silin-like ferromagnetic mode. We discuss the role of the interaction 

1/'^ , parameter Ff on the behavior of the Goldstone mode and the first Silin-like ferromagnetic mode. 

PACS numbers: 71. 10. Ay 

m ■ 

' Landau introduced the Fermi liquid theory in 1956 to describe the low-temperature thermodynamic properties, 

transport and collective phenomena in strongly correlated Fermi systems. The theory of Fermi liquids was used to 
consider spin waves in ferromagnetic metals and further developed into the framework of the ferromagnetic 
. Fermi liquid theory (FFLT). The FFLT can be used to study weak ferromagnetic metals such as MnSi, ZrZn2, 
\^ Ni^Al, subject to increasing interest for their physical properties and potential applications. Many attempts to 
CO i understand these systems focus on the dynamics of the magnetic excitations, a complex problem, both theoretically 
' and experimentally. A direct experimental information on magnetic excitations in ferromagnetic metals can be 
obtained from inelastic neutron scattering measurements. In weak magnets the scattering intensity is relatively low, 
which makes the measurements difficult. The experimental results for Ni^Al are limited |3|, but there are almost 
^ complete measurements for MnSi [Q. 

In this paper we derive the dispersion relations for spin waves in weak ferromagnetic Fermi liquids (at T = 
^ K), using the FFLT, with the following main results: (i) Two modes are obtained - the well known Goldstone mode 
I and a new mode, similar to the first Silin mode in the case of para magnets in a magnetic field which we call 
first Silin-like ferromagnetic mode; (ii) In the hydrodynamic limit (small g-limit, where g = j^fl is the momentum 
transfer) our model leads to the usual quadratic dispersion law uj{q) ^ for the Goldstone mode, but there are 
^ deviations from the hydrodynamic results for the Goldstone mode and for the first Silin-like mode for larger values of 
q; (iii) The Goldstone mode touches the Stoner excitation boundary tangentially, as predicted by the random phase 
approximation (RPA) theory jo), but in contrast with the RPA our model gives different shape for the uj{q) curve in 
' the region close to the particle-hole continuum. 

To study the collective excitations in a ferromagnetic metal we investigate the time evolution of the quasi-particle 
distribution function, given by 2 x 2 matrix in spin space [np{r,t)]af3 = np{r,t)Sa[3 + m.p{r,t)aa[3, where Up is the 
quasi-particle density and mp is the magnetization density. We start with the quasi-classical kinetic equation 

Oti 

-Q^ = [np: ep\p.B. + [rrip, hp]p,B. (1) 

where P.B. stands for Poisson brackets. The evolution of the quasi-particle distribution function is generated by 
the quasi-particle Hamiltonian [ep(f,t)]ai3 = Sp{f,t)5af} + hp(r,t)aai3, where ep is a mean quasi-particle energy and 
hp = —B + 2 '}2ipi Ipp/inpi is an effective magnetic field, which includes both coupling to an external field B and 
to an internal field, due to the quasi-particle interactions (we set the nuclear magnetic moment — 1). The 

interaction parameter /^p, is the antisymmetric component of the quasi-particle interaction fpp,' (in the limit of 
weak ferromagnetism f^p/ can be treated as rotationally invariant in spin space 0). From Eq.(l) we get a set of 
coupled equations, describing the time evolution of the quasi-particle density and the magnetization density. For 
small fluctuations 6np and Snip around the ground state these equations decouple. We are interested only in the 
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equation, which describes oscillations of the transverse components of the total magnetization. The linearized kinetic 
equation for Snip is ^, ^ 



= -2(to° X Shp + 6nip x + I[mp] 



(2) 



where I[mp] is the collision integral, Shp is the fluctuation of the effective magnetic field and hp = J2p' fpp'^'i'p' is the 



equilibrium field due only to the quasi-particle interactions. 

We study the response of the system to a transverse perturbation 6B = SB^i + 5Byj (the transverse spin 
polarization 5rn^ — Sm'p + iSmy, there is an analogous equation for 5mZ^). For excitations close to the Fermi surface 
we introduce the quantities Snip = —{dn^/ de^)v^ and ra^ — —{dnp/dep)mP/N{0), where A^(0) = kpni* /ir'^ is 
the density of states at the Fermi surface and m* is the effective mass of the quasi-particle ^, ^. After Fourier 
transformation of the linearized kinetic equation, Eq.(2), we obtain 



{i^-vp- q + 2m^fS)y+ - (2mO + N{0)vp ■ q) / — /pV^^^ 



dVl' 



f 2mO 



SB^ 



(3) 



We expand both and the Landau parameters / " , in a series of Legendre polynomials. 



For the collision term I[i'p] we use a relaxation-time approximation I[{yp)i\ = 

since quasiparticle spin is conserved 



p „„„ „„„ „„„„„„ ^„»„.„„„„»„ 
and f-p, = Y.j!m-P') 

_1+£l^H^+11(^+)^^ where tu is the spin-diffusion lifetime and I[{Vp)t=Q\ 
in coUisions g]. 

We analyze the possible modes of oscillation using two different approaches: (i) Low-order distortions approxi- 
mation (hydrodynamic model); (ii) All-order distortions model. 

We start with the hydrodynamic model and project out the £ = 0, 1, 2 components of Eq.(3). The Vp ■ q term 



in Eq.(3) couples the distortions of order 
Vp - q^ VFqcos{p ■ q)). 



to distortions of order £ ± 1 (for excitations near the Fermi surface 
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5B^ 



TD 



-u+{l + FS)qvF--i^i 



1 H — — qvp ~ —SB^qvp 



N{0) 
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pa ^2 



iil + Fi/5) 



TD 



2 + 



pa 

i + — jqvF- yi^3 



l + ^]qVF = 



(4) 



(5) 



(6) 



where F° is a dimensionless Landau parameter, defined as F° = N{Q)f^. We take the limit SB^ and for 
simplicity work in the extreme coUisionless regime (T = K). The distortions i/^ can be neglected in the small q- 
limit if qvF <^ \'-^e{q) ^ A|, where i-ui{q) is the dispersion relation for the spin- wave mode considered in any particular 
case and A = —2imPFg/N{Q) is the Stoner gap parameter. From Eq.(4), Eq.(5) and Eq.(6) we obtain the dispersion 
equation 



{UJ - UJ+){UJ - UJ+) 



Fl\ {qvF? 
3 



+ -^(l + ^^o'')(l + f)^ = (7) 



where uj^ — A(l — F^ /'iF^) and = A(l — F2 /hF^) are the eigenfrequencies of the spin-wave modes uj'^{q). In 
what follows we assume F° = for ^ > 1. Then the first solution of Eq.(7) is 



(l + F-)(l + F{'/3) {qVFf 



O {{qvp) 



(8) 



which we identify as the well known quadratic dispersion law for the H. — Q Goldstone mode. In the limit of weak 
ferromagnetism F^ — > — 1~ and the Goldstone mode is almost independent of the interaction parameter F^. The 
second solution of Eq.(7) 



(l + fo°)(l + Ff/3) 4 (l + Ff/3) 
UJ+ 5 I^F'^/iFg 



(qvpf 



(9) 
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gives the dispersion relation for a new spin-wave mode, the i = 1 mode, which we call first Silin-like ferromagnetic 
mode. The second term of order (qvp)^ in Eq.(9) comes in because the distortions 1^2 taken into account; for 
a certain range of values for the interaction parameter this term is comparable in size with the first quadratic 
term, caused by the i^q distortions. The assumption F2 = forces the eigenfrequency of the second Silin-like 
ferromagnetic mode {q) to the Stoner gap (w^ = A) and makes the mode strongly damped. The second Silin-like 
ferromagnetic mode is propagating only if F2 ^ 0. 

Eq.(4) and Eq.(5) are the two hydrodynamic equations for the magnetization density - the net magnetization 
conservation law and the equation of motion for the magnetization current. The equation for the magnetization 
current contains a precession term from the internal field i'^[2'mP /N{0)qvF]F^ /3, found also in the paramagnetic 
case ]lO[ . The precession term in Eq.(5) is missing in the ordinary hydrodynamic equations for ferro magnets 
these equations lead to the quadratic dispersion law for the Goldstonc mode, and miss completely the first Silin-like 
ferromagnetic mode. In addition, the equation for the magnetization current contains distortions of the Fermi surface 
of order £ — 2, z/^. While these distortions produce corrections of higher order in the dispersion relation for the 
Goldstone mode, they are important for the first Silin-like ferromagnetic mode, as shown above. 

At present we do not have data for the value of the interaction parameter F^ in weak ferro magnets, but in other 
strongly correlated Fermi systems like liquid ^He and ^_ffe — He mixtures the range of F^ is |Ff | < 1, We plot 
in Fig.l the dispersion of the Goldstone mode, Eq.(8) and the dispersion of the first Silin-like ferromagnetic mode, 
Eq.(9) for Fq = —1.1 and for different values of F^, |Ff | < 1.5. Changing the value of F^ does not influence the 
behavior of the Goldstone mode, but produces important effects on the first Silin-like ferromagnetic mode. If F° = 0, 
the first Silin-like ferromagnetic mode starts exactly from the Stoner gap. The mode is strongly damped because it 
stays right on the particle- hole continuum, defined d,s A ± qvp. If F^ ^ 0, the first Silin-like ferromagnetic mode 
moves up or down with respect to the Stoner gap, depending on the sign of F^ (see Fig.l), and can propagate before 
entering into the particle-hole continuum and being damped. 

W.^ "I 
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Figure 1: The hydrodynamic modes <^^{q) = 0+ and i^^{q) = 1^ for ferromagnetic Fermi liquid for Fq = — 1.1 and 
= —1.5, —0.5,0, 1.5 (shown in the brackets). 

We make a rough estimate of the coefficient in front of q^ in the dispersion relation for the Goldstone mode, 
Eq.(8), and compare it with the coefficient in front of q^ in the empirical dispersion law for spin- wave excitation for 
MnSi, nujq[meV] = 0.13 + 0.52g2[A-2] pa ^ _2 pa ^ j,^ ^ 3Q0 j^^ ^nd use the data available 



for MnSi: the equilibrium magnetization under a pressure P = 6.6 kbar is mP = 0.389 ^b/Mu |12 and the density 
of states iV(0) = 19 states/Ry Mn atom spin, obtained from band structure calculations jl^. Then Eq.(8) has the 
form LUQ{q) ~ 2.5g^, where Lo^iq) is scaled with respect to the Fermi energy ep and q - with respect to the Fermi 
momentum kp. The empirical law expressed in the same units is uj{q) ~ 2q^, thus we obtain the same order of 
magnitude for both coefficients. 

To analyze the spin-wave modes beyond the small g-limit, we go back to the kinetic equation, Eq.(3), and derive 
the all order distortions model. As before we expand and f^^, in a series of Legendre polynomials and obtain a 
general expression for distortions of order € 
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Pi{x)Pi.{x)dx 



n,. 



PeXx)dx 
N{0)qvF x-s 



(10) 



where x cos{p-q), s = {uj-A)/qvF, f^«'(s) = ^e'e{s) = 5 /_i ^M^lljd^ldx and rioo(s) = 1 + ^sln{s - l)/(s + 1). 
All distortions of order £'= 0, 1, 2, 3,..., 00 enter explicitly Eq.(lO) as long as the interaction parameter ^ 0. 

We are interested in the effects of the Landau parameters Fq and on the dispersion of the Goldstone mode 
and the first Silin-like ferromagnetic mode and as in the hydrodynamic model, we set F° = for ^ > 1. To calculate 
the dynamical response w) of the system to a transverse perturbation SB~^, x(g, cj) = —N{0)vq /SB'^ , we write 



Eq.(lO) for i = and use the net magnetization conservation law. 
The real poles of the response function xio", ^) 



Eq.(4), to express i^i in terms of Vq and SB~^. 
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(11) 



give the undamped spin waves. In the small g-limit of Eq.(ll) we recover the dispersion equation from the hydrody- 
namic model, Eq.(7), with F-f = 0. 

We solve numerically Eq.(ll) for \s\ > 1 and obtain two solutions oJo{q) ^^^d uj^{q), which we identify as the 
Goldstone mode and the first Silin-like ferromagnetic mode. When q increases both modes i^oiq) and uj^{q) start 
to deviate from the hydrodynamic values. Close to the particle-hole continuum (s — > ±1), where the higher order 
distortions become significant, there is a big shift in the frequency (see Fig. 2 and Fig. 3). To see the behavior of the 
modes close to the singularity we expand Eq.(ll) in the small e-limit (s = —1 — e and e 0^) and obtain 



pa 

3 



e{q) 1 + ^ = 1 



pa 



FSil + F^/3) 
1 — A/qvF 



(12) 



qvpF^ 

From Eq.(12) we see that there are two singularities at qo = ~A/F^vf and qi = -qoF^/{l + F^/3 + F^/F^) 
singularity at qo is completely determined by the interaction parameter Fq through the relation mP ~ |1 
and it does not depend on the interaction parameter Ff , while the singularity at qi depends critically on Ff . The 
two singularities qo and qi coincide when Ff = — Fq /(4F^/3 -I- 1) and for qi > qo the first Silin-like ferromagnetic 
mode crosses the Goldstone mode. 



The 

pa^l/2 



0.01 - 



0.00 



0.00 O.Oi 



Figure 2: The mode ujQ'{q) = 0^ with the hydrodynamic mode uj^{q) = 0+ for Fq = —1.1 and Ff = 1.5. 

We shall discuss now the behavior of the Goldstone mode in the vicinity of the singularity at qo . In agreement 
with the RPA ^ the Goldstone mode touches tangentially the particle-hole continuum boundary. But for q > qo 
we find a branching of the solution uio{q) of Eq.(ll), i.e. there are two different frequencies 1^0 i^) ^'^^ ^^'^ same q. 
Close to the singularity at qo, for q ^ qo the asymptotic expansion uj{q) = A — [1 -I- e{q)]qvF describes the behavior 
of the upper brunch of the Goldstone mode W(f (g) (here e ^ 0+, s — > — 1 and Eq.(ll) has real solutions tor q q^ ). 
Thus our model predicts that the Goldstone mode curves upward and then touches tangentially the particle-hole 
continuum. It can be shown that including the higher order interaction parameter Fg does not change the qualitative 
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behavior of the Goldstone mode and the singularity at go (90 does not depend on Ff and F^). Therefore, for qo > qi 
there is no a value of the Stoner gap parameter for which the dispersion curve i^o{q) goes through a maximum, as 
predicted by the RPA 1^. It is worth noting that for values of the interaction parameter < —F^/{4Fq/3 + 1) the 
effect of crossing of the Goldstone and the first Silin-like ferromagnetic modes may result in changing the shape of 
the Goldstone mode in the region of larger momentum transfer, forcing it to curve down and then touch tangentially 
the Stoner excitation boundary, i.e. forcing the dispersion curve uJo {q) to go through a maximum. 
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Figure 3: The mode (jJi {q) = 1^ with the hydrodynamic mode i^iiq) = 1^ for Fq — —1.1 and Ff = —1.5,1.5 
(shown in the brackets) 

When the first Silin-like ferromagnetic mode i^i'{q) approaches the particle-hole continuum, it shows different 
behavior for different values of F^ (see Fig. 3). For negative values of F° the singularity at qi determines the point 
where (g) touches the particle-hole continuum, for Ff = the singularity is at gi = and Wi*(g) degenerates to 
a point, for positive values of F" the mode uj^{q) approaches asymptotically the particle-hole continuum. 

We derived above a new spin-wave mode, the first Silin-like ferromagnetic mode. The mode is propagating for 
any values of the interaction parameter Ff 0. Now the obvious question is why the first Silin-like ferromagnetic 
mode has not been observed by inelastic neutron measurements on weak ferro magnets? To answer to this question, 
we calculate the intensities [ |l^ of the hydrodynamic Goldstone and first Silin-like ferromagnetic modes, extending 
our model to low temperatures and assuming the regime where the spin-diffusion lifetime varies with the temperature 
as TD oc 1/T^. It can be shown that the intensity of the Goldstone mode is about 3 orders of magnitude higher 
than the intensity of the first Silin-like ferromagnetic mode, the linewidth Tq of the Goldstone mode is much less 
than the linewidth Fi of the first Silin-like ferromagnetic mode (the ratio Fi/Fg = 1 + (1 + F^)uj^ / ujq (q) depends 
on the Landau parameters Fq and F^) and with increasing the temperature the intensities of both modes decrease 
and broaden. 

In summary: our main results, the dispersion laws for the Goldstone mode and for the first Silin-like ferromagnetic 
mode, are obtained as an implication of the FFLT. Including the interaction parameter F° brings in new physics 
- the presence of a new spin wave mode, the first Silin-like ferromagnetic mode. The framework presented in this 
paper allows us to obtain the whole picture for the dispersion and attenuation of the spin waves in ferromagnetic 
metals. 
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